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Abstract . This report is the first of a series dealing with the 

determination of optimal trajectories for the aeroassisted flight experiment 

J* J 

(AFE) I The AFE refers to the study of the free flight of an autonomous 
spacecraft, shuttle-launched and shuttle-recovered. Its purpose is to gather 
atmospheric entry environmental data for use in designing aeroassisted 
orbital transfer vehicles (AOTV). 

It is assumed that: 'ihe spacecraft is a particle of constant mass;: ) 

the Earth is rotating with constant angular velocity; the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored; 
the atmosphere is at rest with respect to the Earth. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also 
for AOT problems and space shuttle problems) are derived in an Earth-fixed 
system. Transformation relations are supplied which allow one to pass from 

quantities computed in an Earth-fixed system to quantities computed in an 

, ■ \ 

inertial system, and viceversa. 5 

Key Words . Flight mechanics, hypervelocity flight, atmospheric flight, 
coordinate systems, equations of motion, transformation techniques, optimal 
trajectories, aeroassisted flight experiment, aeroassisted orbital transfer, 
space shuttle reentry. 
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1 . Introduction 

This report is the first of a series dealing with the determination 
of optimal trajectories for the aeroassisted flight experiment (AFE). The 
AFE refers to the study of the free flight of an autonomous spacecraft, 
shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric 
entry environmental data for use in designing aeroassisted orbital transfer 
vehicles (AOTV). 

It is assumed that: (a) the spacecraft is a particle of constant 
mass; (b) the Earth is rotating with constant angular velocity; (c) the 
atmosphere is at rest with respect to the Earth; (d) the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also 
for AOT problems and space shuttle problems) are derived in an Earth-fixed 
system. Transformation relations are supplied which allow one to pass from 
quantities computed in an Earth-fixed system to quantities computed in an 
inertial system, and viceversa. 

Previous Research . Previous research on the topics covered here can be 
found in Refs. 1-11. For the general theory of flight paths and coordinate 
systems, see Refs. 1-2; for the equations of flight over a spherical Earth, 
see Refs. 1-3; for the perturbed motion about an oblate Earth, see Ref. 4; 
for AFE problems, see Ref. 5; for reentry problems, see Ref. 6; for methods 
of orbit determination, see Refs. 7-8; for the values of the astrophysical 
quantities, see Ref. 9; for the values of the characteristic constants of 
the oblate Earth, see Refs. 10-11. 
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Outline . Section 2 contains the notations, and Section 3 defines the 
basic coordinate systems. The relations between coordinate systems are 
discussed in Section 4, and the angular velocity (or evolutory velocity) 
is introduced in Section 5. The kinematical equations for an Earth-fixed 
system are derived in Section 6, and the dynamical equations are obtained 
in Section 7. Section 8 summarizes the results, and Section 9 presents the 
transformation relations which allow one to pass from quantities computed in 
an Earth-fixed system to quantities computed in an inertial system, and viceversa. 


.laKlilub. iJ . Mi i i|tP § n.Mi lii 
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2. Notations 

Throughout the paper, the following notations are employed: 
a = acceleration, m/sec 2 ; 

A = aerodynamic force, N; 

Cp = drag coefficient; 

C L = lift coefficient; 

Cq = side force coefficient; 

D = drag force, N; 

2 

f = latitudinal component of the gravitational acceleration, m/sec ; 

2 

g = radial component of the gravitational acceleration, m/sec ; 

J 2 = characteristic constant of the Earth's gravitational field; 

J 3 = characteristic constant of the Earth's gravitational field; 

= characteristic constant of the Earth's gravitational field; 

L = lift force, N; 

m = mass, kg; 

M = Mach number; 

Q = side force, N; 

r = radial distance, m; 

r g = equatorial radius, m; 

r = polar radius, m; 

P 

R g = Reynolds number; 

2 

S = reference surface area, m ; 

T = thrust force, N; 

2 2 

U = Earth's gravitational potential, m /sec ; 

V = velocity, m/sec; 
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W = gravitational force, N; 

x = Cartesian coordinate, m; 
y = Cartesian coordinate, m; 
z = Cartesian coordinate, m; 
a = angle of attack, rad; 
y = path inclination, rad; 

0 = longitude, rad; 
y = bank angle, rad; 

3 2 

y g = Earth s gravitational constant, m /sec ; 

3 

p = air density, kg/m ; 
a = sideslip angle, rad; 

<J> = latitude, rad; 

X = heading angle, rad; 

w = angular velocity of the Earth with respect to an inertial system, rad/sec 

* angular velocity of the local horizon system with respect to the Earth 

axes system, rad/sec. 

Subscripts 

b = body axes system; 
e = Earth axes system; 
h = local horizon system; 

1 = inertial system; 

w = wind axes system. 

Superscripts 


derivative with respect to time; 
vector quantity. 
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3. Basic Coordinate Systems 

The basic coordinate systems for flight over a spherical Earth are the Earth 
axes system Ox g y e z e , the local horizon system Px b y b z b > the wind axes 
system Px^y^, and the body axes system Px b y b z b . 

3.1. Earth Axes System . The Earth axes system Ox g y e z e is a Cartesian 
reference frame which is rigidly attached to the Earth. Its origin 0 is 
the center of the Earth; the z g -axis is aligned with the axis of rotation of 
the Earth and is positive northward; the axes * e ,y e are orthogonal to the 

z -axis and are directed radially; the trihedral Ox y z is right-handed. 

0 c c e 

In particular, the plane x g ,y e contains the fundamental parallel (the Equator); 
and the plane x g ,z e contains the fundamental meridian (the Greenwich meridian). 
The symbols ^ e >J e **< e denote the unit vectors of the Earth axes system. 

3.2. Local Horizon System . The local horizon system Px h y h z h is a Cartesian 
reference frame defined as follows. Its origin P is identical with the 
instantaneous position of the spacecraft; the z h ~axis is directed radially 
(that is, vertical) and is positive downward; the axes x h ,y h are orthogonal 

to the z h -axis (therefore, they are tangent to the spherical surface through 
P; they form the so-called local horizon plane); the trihedral P x b y b z ^ is 
right-handed. In particular, the x. -axis is tangent to the local parallel 
through P and is positive eastward; the y h -axis is tangent to the local 

-v ' + 

meridian through P and is positive southward. The symbols i h , ,k h denote 
the unit vectors of the local horizon system. 

3.3. Wind Axes System . The wind axes system P x w y w z w is a Cartesian 
reference frame defined as follows. Its origin P is identical with the 
instantaneous position of the spacecraft; the x^-axis is tangent to the 
flight path (relative velocity) and is positive forward; the axes y w » z w 
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are orthogonal to the x w ~axis and are such that the trihedral 
is right-handed. In particular, the z -axis is contained in the plane of 
symmetry of the spacecraft and is positive downward for the normal flight 
attitude of the spacecraft; the y w -axis is positive rightward for the 
normal flight attitude of the spacecraft. The symbols i w >J w > k w denote 
the unit vectors of the wind axes system. 

3.4. Body Axes System . The body axes system P^y^ is a Cartesian 
reference frame defined as follows. Its origin P is identical with 
the instantaneous position of the spacecraft; the y b ~axis is orthogonal to 
the plane of symmetry of the spacecraft and Is positive rightward; the 
axes x b » 2 b are orthogonal to the y^-axis, are contained in the plane of 
symmetry, and are such that the trihedral is right-handed. In 

particular, the x b ~axis is positive forward, the y b -axis is positive 
rightward, and the z b -axis is positive downward for the normal flight 
attitude of the spacecraft. The symbols i b ® J b’^b denote the unl 't vec tors 


of the body axes system. 
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4. Relations between Coordinate Systems 

In this section, the relationships between the different coordinate 
systems are derived; more specifically, attention is focused on the 
following system pairs: Earth axes-local horizon; local horizon-wind axes; 
and wind axes-body axes. 

We recall that, in the Earth axes system, a point P can be described 
via its Cartesian coordinates x e ,y e ,z g . Alternatively, P can be described 
via its spherical coordinate r, 6, <J>. Here, r is the radial distance from 
the center of the Earth; 6 is the longitude, positive eastward; and <f> is 
the latitude, positive northward. 

4.1. Transformation from Earth Axes to Local Horizon . The local 
horizon system Px^y^z^ can be obtained from the Earth axes system Ox e y g z e 
by means of the combination of four rotations and one translation. This 
requires the definition of four intermediate coordinate systems: the system 
0 x ^1 z-j ; the system 0x 2 y 2 z 2 ; the system Px-^z^, and the system Px 4 y 4 z 4 » 

The system Ox-^z-j is obtained from the Earth axes system 0x g y e z e by 
means of the counterclockwise rotation 6 around the z e ~axis. Note that the 
z-j-axis is the same as the z e -axis, that the axes x-j ,y-j are contained in the 
equatorial plane, and that the axes x ] ,z 1 are contained in a meridian plane. 
The symbols i, ,k 1 denote the unit vectors of the system Ox-jy-j z-j . 

The system 0x 2 y 2 z 2 is obtained from the system Ox^-jZ^ by means of 
the clockwise rotation 4> around the y-j-axis. Note that the y 2 ~axis is 
the same as the y-j-axis, that the axes x 2 ,z 2 are contained in a meridian 
plane, and that the axes y 2 ,z 2 are contained in a plane parallel to the local 
horizon plane. The symbols i 2 ,j 2 »k 2 denote the unit vectors of the system 
Ox 2 y 2 z 2 • 
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The system Px 3 y 3 z 3 is obtained from the system Ox 2 y 2 z 2 by means of 
the radial translation r, leading from point 0 to point P. Since there 
is no rotation, the axes x 3 ,y 3 ,z 3 are parallel to the axes x 2 ,y 2 ,z 2 ; 
in particular, the axes x 3> z 3 are contained in a meridian plane, while 
the axes y 3 ,z 3 are contained in the local horizon plane. The symbols 
denote the unit vectors of the system Px 3 y 3 z 3 . 


The system Px^y^z^ is obtained from the system Px 3 y 3 z 3 by means of 

the counterclockwise rotation tt/ 2 around the z 3 ~axis. Note that the 

z^-axis is the same as the z 3 -axis, that the axes y 4 » z 4 are contained in a 

meridian plane, while the axes x 4 ,z 4 are contained in the local horizon plane. 

-v -y 

The symbols i 4 » J 4 * k 4 denote the unit vectors of the system Px 4 y 4 z 4< 

The local horizon system Px^z^ is obtained from the system Px 4 y^z 4 by means 


of the clockwise rotation tt/2 around the x 4 ~axis. Note that the x^-axis 
is the same as the x 4 ~axis, that the axes y h * z h are contained in a meridian 
plane, and the axes x^y^ are contained in the local horizon plane. 

In vector-matrix notation, the successive transformations leading from 
one coordinate system to another can be expressed as follows: 



COS0 

sin0 

0 


-4- 

-sinB 

COS0 

0 


j e 

0 

0 

1 


-> 

k 



• 


e 

— 




■ 

cos<f> 

0 

s i n<p 


j i 

0 

1 

0 


-> 

j i 

-sincj) 

0 

cos<j> 


M 


(la) 


(lb) 
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while Eqs. (2) imply that 



(2a) 


(2b) 


(2c) 


(3) 


Therefore, upon combining Eqs. (3)-(4), we see that 


* 
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T 


-sine 


cose 

0 

-y 

= 

cose 

s i n4> 

sine sinc|> 

-COS(f> 

sz 
t ^ 

1 


-cose 

COS<j) 

-sine cos4i 

-sin<f> 


(5a) 


with the implication that 


1 e 


-sine 

cose sin<j> 

-cose cos<J> 

->► 

j e 

= 

cose 

sine sincf> 

-sine cos<|> 

-y 

k e J 


0 

-COS<t> 

-si nd> 

4.2. 

Transformation 

from Local Horizon to Wind 



■ -y 




-y 


j h 


-y- 


k h j 


(5b) 


system P x v/ y w z w can be obtained from the local horizon system Px^y^h by 
means of the combination of three rotations. This requires the definition 
of two intermediate coordinate systems: the system PXgygZg and the system 

Px 6 y 6 z 6‘ 

The system Px 5 y 5 z 5 is obtained from the local horizon system Px h y h z h 
by means of the counterclockwise rotation x around the z h ~axis. Note 
that the z^-axis is the same as the z^-axis, that the axes Xg,yg are contained 
in the local horizon plane, and that the axes Xg,Zg are contained in the 
plane (OP ,V), where OP is the radius vector connecting the points 0 and P 
and V is the spacecraft velocity vector. Also note that the axis Xg has the 

•f 

direction of the projected velocity vector V p ; this is the projection of 
V on the local horizon. The angle x 1S called the heading angle and is 

—y ^ 

positive if the projected velocity vector V p is directed outward with respect 
to the local parallel. The symbols ig,jg,kg denote the unit vectors of the 


system PXgygZg. 
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v 


H. 


The system Px^y^z^ is obtained from the system Px g ygZg by means of 
the counterclockwise rotation y around the yg-axis. Note that the y g -axis 
is the same as the yg-axis and that the axes x g ,z g are contained in the 

*>■ ~y 

plane (0P,V). Also note that the x g -axis is positive forward and that the 
Zg-axis is positive downward. The angle y is called the path inclination 

-y 

and is positive if the velocity vector V is inclined upward with respect 

-y -y ~y 

to the local horizon. The symbols ig^gjkg denote the unit vectors of the 
system PxgygZg. 

The wind axes system Px^^ is obtained from the system PXgy g Zg 
by means of the counterclockwise rotation p around the Xg-axis. Note that 
the x -axis is the same as the x c -axis. Also, note that the x-axis is 
positive forward, the y w ~axis is positive rightward, and the z w -axis is 
positive downward and is contained in the plane of symmetry of the spacecraft. 
The angle p is called the angle of bank and is positive if the spacecraft 
is banked to the right. 

In vector-matrix notation, the successive transformations leading 
from one coordinate system to another can be expressed as follows: 



cosy 

siny 

* 

0 


■ - 

-y 

siny 

cosx 

0 


j h 

-y 

0 

0 

1 


> 





r -y - 

cosy 

0 

-siny 


i 5 

-y 

0 

1 

0 


h 

-y 

siny 

0 

cosy 


k 5 


(6a) 


(6b) 
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1 

0 

0 


0 

cosy 

-siny 


■ 


■ 

0 






siny 






cosy 


k 6 

. 


. J. 


(6c) 


Equations (6) lead to 


T ->- 

i 

w 


cosy cosy C0S Y sinx -siny 


= 

siny siny cosy siny siny sinx siny cosy 
-cosy siny +cosy cosx 

K 


cosy siny cosy cosy siny sinx cosy cosy 
+siny siny -siny cosx 



(7a) 


with the implication that 


-V 

iu 


cosy cosy 
cosy siny 
-siny 


siny siny cosx 
-cosy siny 

siny siny sinx 
+cosy cosx 

siny cosy 


cosy siny cosy 
+siny siny 

cosy siny sinx 
-siny cosx 

cosy cosy 



. (7b) 


4.3. Transformation from Wind Axes to Body Axes . The body system 

Px b y b z b can be obtained from the w1nd axes s y stem Px w y w Z w by means the 
combination of two rotations. This requires the definition of one intermediate 

coordinate system, the system Px 7 y 7 z 7 . 

The system Px 7 y ? z 7 is obtained from the wind axes system by 

means of the counterclockwise rotation a around the z w ~axis. Note that the 

z 7 -axis is the same as the z w -axis and that the axes x y ,z 7 are contained 

in the plane of symmetry of the spacecraft. Also note that the axis x 7 
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is positive forward, the axis y 7 is positive rightward, and the axis z 7 
is positive downward. The angle a is called the sideslip angle and is 
positive if the velocity vector V is directed leftward with respect to the 
plane of symmetry of the spacecraft. 

The body axes system Px^y^z^ is obtained from the system Px 7 y 7 z 7 by means of 
the counterclockwise rotation a around the y 7 ~axis. Note that the y^-axis 
is the same as the y 7 -axis and that the axes x b ,z b are contained in the 
plane of symmetry of the spacecraft. Also note that the axis x b is positive 
forward, the axis y b is positive rightward, and the axis z b is positive 
downward. The angle a is called the angle of attack and is positive 
if the velocity vector V is directed downward with respect to the x b ~axis 


of the spacecraft. 

In vector-matrix notation, the successive transformations leading 
from one coordinate system to another can be expressed as follows: 


■ -> 


* 


■ ' 

V 


cosa 

sina 

0 

-+ 




0 

h 

= 

-sina 

cosa 



0 

0 

1 



- 


- 

_ . 


_ 










cosa 

0 

-sina 

^b 

= 

0 

1 

0 



sina 

0 

cosa 

I m m 


. 


1 


-V 



with the implication that 
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T 





Jb 

II 



k b 



cosa 

cosa 

cosa 

sina 

-sina 

-sina 


cosa 


0 

sina 

cosa 

sina 

sina 

cosa 


T 

" 


'w 


->■ 


j w 


-+ 


k 

w 

.1 

■ « ► 


and that 



cosa cosa 
cosa sina 
-sina 


sina 

sina 

cosa 

cosa 

sina 

sina 

0 

cosa 




(9a) 


(9b) 
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5. Angular Velocity 

In this section, we compute the angular velocity (or evolutory velocity) 
of the local horizon system with respect to the Earth axes system. To do 
so, consider the behavior of the spacecraft between the time instants t and 

- 4 - 

t + dt, and denote by d^ he the infinitesimal vectorial rotation of the local 
horizon system with respect to the Earth axes system. This infinitesimal 
vectorial rotation can be decomposed into partial rotations as follows: 


d ^he 

= dfi^ + dQ^^ + <^32 + d^i + dfi-j ^ , 

(10) 

d ^h4 

= o. 

(Ha) 

dn 43 

= o. 

(lib) 

dfi 32 

= 0, 

(He) 

d^i 

=-d<J> j-| > 

(lid) 


II 

CL 

CD 

7n- ± 
CO 

• 

(He) 


Here, de denotes the infinitesimal change of the longitude and d<{> denotes 
the infinitesimal change of the latitude. Note that the rotation d0 occurs 
around the z g -axis and is positive counterclockwise and that the rotation 
d4> occurs around the y-j-axis and is positive clockwise. This explains the 
difference in the signs appearing on the right-hand sides of Eqs. (lid) 
and (lie). 
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Upon combining Eqs. (10)-(11), we see that the infinitesimal vectorial 
rotation of the local horizon system with respect to the Earth axes system 
can be written as 


d^ e = d9 k g - dtj> j’i . (12) 

As a consequence, the angular velocity of the local horizon system with 
respect to the Earth axes system is given by 

“he = d ^he /dt = ~ (13) 

In the light of Eqs. (l)-(5), the unit vectors k g and j-j can be expressed 
in terms of the unit vectors of the local horizon system as follows: 


k g = -cos<{> j h - sinc}> k h , 
J 1 = 1 h’ 


(14a) 

(14b) 


so that 


“he = ”^h " ® cos<t> ^h ‘ ® sin<J) k h’ 

Next, Poisson's formulas are employed to compute the derivatives of 
the unit vectors of the local horizon system with respect to time: 


di^/dt = co 

dj h /dt = co 
->■ 

dk^/dt = a) 


he 

he 

he 





(16a) 

(16b) 


(16c) 
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Upon combining Eqs. (15)-(16), we obtain the relations 
di h /dt = - (esin<J>)j h + (0cos<l>)k h , 
dj h /dt = (0sin<|>)i h - 4>k h » 

dk^/dt = - (ecos<}>)i h + 4>J h , 


(17a) 

(17b) 

(17c) 


whose vector-matrix form is the following: 



-y 

-y 


• • 

0 -0sin4> 6cos4> 

• « 


-y 

n 'h 

-y 

(d/dt) 

j h 

.y 


0sin<J> 0 -<f> 

• • 

-9cos({) 4> 0 

. * 


j h 

.y. 


(18) 


It is interesting to note that Eq. (18) can also be obtained by taking 
the time derivative of Eq. (5a) and using Eq. (5b). 
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6. Kinematical Equations 

In this section, we derive the scalar relationships corresponding to 
the vectorial equation 

dOP/dt = V. (19) 

-V 

Here, V denotes the velocity of the spacecraft with respect to the Earth 
and OP denotes the position vector joining the center of the Earth 0 with 
the spacecraft position P. 

First, we observe that the position vector OP is given by 

OP * ~rk h , (20) 

-F 

where r is the radial distance from the center of the Earth and is 
the third unit vector of the local horizon system. As a consequence, the 
time derivative of OP can be written as 

dOP/dt = -rk h - r(dk h /dt), (21) 

where dk^/dt is given by Eq. (17c). Therefore, upon combining Eqs. (17c) 
and (21), we obtain the relation 

dOP/dt = 0rcos<|) 1. - $rj h - rk^. (22) 

Next, we observe that the velocity vector V is given by 

V = vT , (23) 

w 

where V is the velocity modulus and i is the first unit vector of the 

w 

— V 

wind axes system. In the light of Eq. (7a), the unit vector i w can be 


written as 
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-V -* 

i = cosy cosx i h + cosy siny J h - siny k h> (24) 

Therefore, upon combining Eqs.(23) and (24), we obtain the relation 

-> + “► . . 

V = Vcosy cosx 1 h + V cosy siny j' h - Vsiny k h . (25) 

Finally, upon combining Eqs. (19) , (22), and (25), and upon projecting 
the resulting vectorial equation on the axes of the local horizon system, 
we obtain the following scalar form of the kinematical equations: 


• 

6 = Vcosy cosy/rcosef). 

(26a) 

• 

<J> =-V cosy siny/r. 

(26b) 


r = Vsiny. 


(26c) 
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7. Dynamical Equations 

In this section, we derive the scalar relationships corresponding 
to the vectorial equation 

T + A + W = ma .j , (27) 

where T is the thrust, A is the aerodynamic force, W is the gravitational 

-V 

force, m is the mass of the spacecraft, and a^ is the inertial acceleration. 

We consider the case where the engine is shut-off, so that 

T = 0, (28) 

and the mass of the spacecraft is constant. Hence, Eq. (27) is written as 

A + W = ma.j. (29) 

Because of the theorem of composition of the accelerations, the inertial 
acceleration can be written as the sum of the relative acceleration dV/dt 

(acceleration with respect to the Earth), the Coriolis acceleration 2w x V, 

->• ^ 

and the transport acceleration oj x (w x OP): 

a. = dV/dt + 2u> x V + Z x (ui x OP). (30) 

Here, to is the angular velocity of the Earth with respect to an inertial 
system. Note that u> is constant and is aligned with the axis of rotation 
of the Earth. Therefore, upon combining Eqs. (29)-(30), we obtain the 
vectorial equation 

A + W = m[dV/dt + 2io x V + u x (m x OP)] . 


(31) 
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We now compute the components of the vectors appearing in Eq. (31) on the axes of 
the local horizon system. 

7,1. Aerodynamic Force . The components of the aerodynamic force on the wind 
axes are the drag D, the side force Q, and the lift L. Therefore, the 
aerodynamic force can be written as 





(32) 


No special significance is implied in the signs appearing on the right-hand 
side of Eq. (32). These signs merely reflect the conventions adopted in 
this report with regard to the positive values for the drag, the side force, 
and the lift. 

In vector-matrix form, Eq. (32) can be rewritten as follows: 



D 0 
0 Q 
0 0 


0 

0 

L 



(33) 


Therefore, upon combining Eqs. 


(7a) and (33), we obtain the relation 


-Dcosy cosy 

-Dcosy siny 

Dsiny 


i 

]r + 

1 

-Qsiny siny cosy 

-Qsiny siny siny 

-Qsiny cosy 


j h 

+Qcosy siny 

-Qcosy cosy 




-Lcosp siny cosy 

-Lcosy siny siny 

-Lcosy cosy 


k h 

-Lsiny siny 

+Lsiny cosy 
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7.2. Gravitational Force . Here, we assume that the Earth is an oblate planet 
and that its mass has radial symmetry with respect to the axis of rotation. 

Because the equatorial radius r is larger than the polar radius r , the 
gravity force W has two components: the radial component mg, directed toward 
the center of the Earth, and the latitudinal component mf, tangent to the 
local meridian and directed toward the Equator. Therefore, the gravity 
force can be written as 

W = mfj h + mgk h . (35) 

The radial component g and the latitudinal component f of the 
acceleration of gravity are related to the Earth's gravitational potential 


U by the expressions 

g = 3U/9r , f = (l/r)3U/3<t>, (36) 

where 

U =-(u e /r)[HJ 2 (r e /r) 2 H 2 + ^(^/r)^ + J 4 (r e /r) 4 H 4 ] , (37a) 

H 2 = 1/2 -(3/2)sin 2 (f), (37b) 

H 3 = (3/2)sincp - (5/2)sin 3 4>, (37c) 

H 4 =-(3/8) + (30/8)sin 2 4> - ( 35/8 ) si n 4 <|>. (37d) 


Here.y is the Earth's gravitational constant, r is the equatorial radius, 
e e 

and J 2 , J 3 , J 4 denote the characteristic constants of the Earth's 
gravitational field. Note that the expression for U is approximate, since 
harmonics of order higher than four are ignored. 
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Upon combining Eqs. ( 36 ) - ( 37 ) , we see that the components of the 


acceleration of gravity can be written as 

g = (y e /r 2 )[l +J 2 (r e /r) 2 6 2 + J 3 (r e /r) 3 G 3 + J 4 (r e /r) 4 G 4 l , (38a) 

G 2 = 3/2 - (9/2)sin 2 cf>, (38b) 

G 3 = 6sin<{> - 10sin 3 <t>, (38c) 

G 4 = -15/8 + (150/8)sin 2 <{> - ( 175/8) sin 4 <f>, (38d) 

and 

f =(u e /r 2 ) [J 2 (r e /r) 2 F 2 + J 3 (r e /r) 3 F 3 + J 4 (r e /^) 4 F 4 ] , (39a) 

F 2 = 3s i n<f) cos<{>, (39b) 

F 3 =-(3/2)coscj) + (15/2)sin 2 <j> cosep , (39c) 

F 4 =-( 1 5/2 ) si n<p cos4> + (35/2)sin 3 <j) cos<|). (39d) 


7.3. Relative Acceleration . Let V xh , V yh , V zh denote the components of 
the relative velocity on the local horizon system. 


V xh = Vcosy cosx, 

V yh = Vcosy siny, 

V zh =-Vsiny. 

With this understanding, the relative velocity (25) can be rewritten as 


° Vh 


V yhY 


V zh k h' 


(40a) 

(40b) 

(40c) 


(41) 
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Therefore, the relative acceleration is given by 

dJ/dt - Vh + Vh + 

+ V xh (d V dt) + V d V dt > + V zh< d V dt >- 

If we combine Eqs. (17) and (26), the time derivatives of the unit 
vectors of the local horizon system can be written as 

di h /dt = -(Vcosy cosy tar»4>/r) j* h + (Vcosy cosx/r)k h , 

-y 

dj h /dt = (Vcosy cosy tan<[>/r)i h + (Vcosy siny/r)k h , 
dk h /dt = -(Vcosy cosy/r)i h - (Vcosy siny/rjj^. 

Upon combining Eqs. (40), (42), (43), the relative acceleration becomes 

dV/dt = (V xh + V 2 cos 2 y cosy siny tan<J>/r + V 2 cosy siny cosy/r)i h 

,2 2 


(42) 


(43a) 

(43b) 

(43c) 


O O 

+ (V . - V^cos^y cos y tancf>/r + V cosy siny si'ny/r)j^ 
+ (V zh + V 2 cos 2 y/r)k h . 


(44) 


7.4. Coriolis Acceleration . The angular velocity of the Earth with respect 
to an inertial system is given by 


to = wk . 
e 


(45) 


Here, k is the third unit vector of the Earth axes system, which is given 
e 

by [see Eq. (5b)] 


k = -cos<f> j. - sin<j> k h . 


(46) 
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Hence, Eq. (45) becomes 

a) = -wcos<}) j^ - usintf) k^. 


(47) 


Next, we recall Eq. (25), 

-> ^ “► . . 

V = Vcosy cosy i^ + Vcosy siny - Vsiny k^. (48) 

Therefore, the Coriolis acceleration becomes 

2w x V = 2coV(siny cos4> + cosy siny sin^>)i h 

- 2ojVcosy cosy si n<J> + 2wVcosy cosy cos4> k^. (49) 

7.5. Transport Acceleration . We recall that the vector connecting the center 
of the Earth with the instantaneous position of the spacecraft is given by 
[see Eq. (20)] 


OP = -rk h . (50) 

We also recall that the angular velocity of the Earth is given by [see Eq. (47)] 

w = -o)cos<j> j. - cosing k^. (51) 

Therefore, the transport acceleration becomes 

oj x ( w x OP) = -w r cos<J) si n<p + to rcos <j> k^. (52) 

7.6. Scalar Equations . Next, we combine Eqs. (31), (34), (35), (44), (49), 
(52). Upon projecting the resulting vectorial equation on the axes of the 
local horizon system, we obtain the following scalar equations: 
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V xh = "( D / m ) C0S 'Y C0S X + (Q/m)(cosy sinx - siny siny cosy) 

- (L/m)(siny siny + cosy siny cosy) 

n n 

- ( V /r ) ( cos y cosy siny tan<j> + cosy siny cosy) 


- 2ojV ( siny cos<}> + cosy siny si n4>) , 

(53a) 

Vyh = -(D/m)cosy siny - (Q/m)(cosy cosy + siny siny siny) 


+ (L/m)(siny cosy - cosy siny siny) + f 


+ (v /r)(cos y cos x tan<j> - cosy siny siny) 


2 

+ 2a)Vcosy cosy sm<f> + to rcos<j> sin<J>, 

(53b) 

V zh = (D/m)siny - (Q/m)siny cosy - (L/m)cosy cosy + g 


9 9 2 2 

- (V /r)cos y - 2wVcosy cosy cos<|> - to rcos <)>. 

(53c) 

We recall that the components of the relative velocity on the axes 
local horizon system are given by [see Eqs. (40)] 

of the 

V xh = Vcosy cosy. 

(54a) 

V yh = Vcosy siny, 

(54b) 

v z h =-Vsiny, 

(54c) 


with the implication that 
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5 xh' 


cosy cosy -siny cosy -siny 


' • " ’ 

V 

V 

= 

cosy siny -siny siny cosy 


Vy 



-siny -cosy 0 


Vcosy y 

m 


w m 


. 


and that 




- 


- • 

V 


cosy cosy cosy sinx -siny 


V xh 

Vy 

= 

-siny cosy -siny siny -cosy 


V 

Vcosy y 


-siny cosy 0 



J 


. ■ 


■ 


( 55 ) 


(56) 


The final step consists of combining Eqs. (53) and (56). This leads 
to the following scalar form of the dynamical equations: 


V = -D/m - gsiny + fcosy siny 

+ io*V(siny cos% + cosy siny cos<j> si n<p) , (57a) 

Vy = (L/m)cosy + (Q/m)siny + (V /r - g)cosy - fsiny siny 

r\ ry 

+ 2a)Vcosy cos4> + to r(cosy cosS - siny siny cost)) si n<p ) , (57b) 

® 9 9 

Vcosy y = (L/m)siny - (Q/m)cosy + (V /r)cos y cosy tan4> + fcosy 

2 

+ 2o)V(cosy sin<)> + siny siny cos<|)) + to rcosy cos<)) sin4>, (57c) 


which can be rewritten as 


V =-D/m - gsiny + fcosy siny 


+ co r(siny cos^cp + cosy siny cost)) sin<t>), 


(58a) 
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y = (L/mV)cosy + (Q/mV)sinu + (V/r - g/V)cosy - (f/V)s1ny sin x 

+ 2wcosx cosef. + (w 2 r/V)( cosy cos 2 c(> - siny sin x cos<}> sin<}>), 

x = (L/mV)siny/cosy - (Q/mV)cosy/cosy 

+ (V/r)cosy cos x tan4> + (f/V)cos x /cosy 

2 

+ 2to(sin(f) + tany sin x cos<t>) + (o> r/V)cos x cos<t> sin<t>/cosy. 


(58b) 


(58c) 
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8. Summary of Results 

In this report, we have derived the equations of motion of a spacecraft 
under the following assumptions: (a) the spacecraft is a particle of constant 
mass; (b) the Earth is rotating with constant angular velocity; (c) the 
atmosphere is at rest with respect to the Earth; (d) the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored. 

An Earth-fixed system has been used, and the following kinematical 
and dynamical equations have been obtained: 


• 

6 = V cosy cosy/rcosej). 

(59a) 

• 

4> =-Vcosy siny/r, 

(59b) 

• 

r = Vsiny, 

(59c) 


and 

V = -D/m - gsiny + fcosy siny 

9 9 

+ go r(siny cos^4> + cosy siny cos<j> sin<f>), (60a) 

y = (L/mV)cosy + (Q/mV)siny + (V/r - g/V)cosy - (f/V)siny siny 

9 9 

+ 2wcosy cos<|> + (w r/V)(cosy cos <j> - siny siny cos<}> sin<J>), (60b) 

y = (L/mV)siny/cosy - (Q/mV) cosy /cosy 

+ (V/r)cosy cosy tan<f> + (f/V)cosy/cosy 

2 

+ 2u)(sin4> + tany siny cos<f>) + (to r/V) cosy cos<J> sin<J>/cosy. 


(60c) 


» Mi i tuiii 
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8.1. Aerodynamic Force. In Eqs. (60), the drag, the side force, and the 
lift are given by 


*** 

D = (l/2)C D pSV 2 , 

(61a) 

— 

Q = (1/2)Cq P SV 2 , 

(61b) 


L = (1/2)C l pSV 2 , 

(61c) 


where Cp is the drag coefficient, Cq is the side force coefficient, 


is the lift coefficient, p is the air density, and S is a reference surface 
area. In turn, the aerodynamic coefficients are functions of the form 


— 

C D = C D^ a5a,M,R e^ 

(62a) 


Cq = Cq(a ,a,M, R q ^ 5 

(62b) 


C L = C^(a,a, M,R g ) , 

(62c) 


where a is the angle of attack, a is the sideslip angle, M is the Mach 
number, and R g is the Reynolds number. 

8.2. Gravitational Force. In Eqs. (60), the radial component and the 
latitudinal component of the acceleration of gravity are given by 

(63a) 
(63b) 
(63c) 


g = (u e /r 2 )H +J 2 (r e /r) 2 G 2 + J 3 (r e /r) 3 G 3 + J 4 (r e /r) 4 G 4 l , 
G 2 = 3/2 - (9/2)sin 2 cf>, 

3 

G 3 = 6sin<{> - lOsin <j>, 

G 4 =— 1 5/8 + (150/8)sin 2 cJ> - (175/8)sinV 


(63d) 
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and 

f =(u e /r 2 )[J 2 (r e /r) 2 F 2 + J 3 (r e /r) 3 F 3 + J 4 (r e /r)\] , (64a) 

F 2 = 3sin<f> cos<J>, (64b) 

=-(3/2)cos<j) + (15/2)sin 2 <}) cost)), (64c) 

F^ =- ( 1 5/2 ) si ncf> costj> + (35/2)sin 3 <J> cos<f>. (64d) 

8.3. Physical Constants . The major physical constants appearing in 
the system (59) -(64) have the following values: 

to = 0.729211595 E-04 rad/sec, (65a) 

y e = 0.39860064 E+15 m 3 /sec 2 , (65b) 

J 2 = 0.10826271 E-02, (65c) 

J 3 =-0.25358868 E-05, (65d) 

J 4 =-0.1624618 E-05, (65e) 

r = 0.6378164 E+07 m, (65f) 

e 

r = 0.6356755 E+07 m, (65g) 

P 

Here, to is the Earth's angular velocity;y is the Earth's gravitational 

e 

constant; J 2 » J 3 > J 4 are the characteristic constants of the Earth's gravitational 

field; r is the Earth's equatorial radius; and r is the Earth's polar radius, 
e P 

Note that the Earthy's sea-level radius r $£ varies with the latitude 4> according 


to the relation 
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r « * < 1/2)(r e + 

V + 

(l/2)(r e - r p )cos(2(J)). 

(66) 

8.4. Spacecraft Data. 

For the AFE vehicle, it is assumed that 


m = 0.16782918 

E+04 

kg. 

(67a) 

S = 0.14314 

E+02 

m 2 . 

(67b) 

a = 0.17000 

E+02 

deg. 

(67c) 

C L =-0.370696 

E+00, 


(67d) 

C D = 0.131452 

E+01 . 


( 67 e ) 


Here, m is the spacecraft mass at atmospheric entry; S is the reference 
surface area; a is the angle of attack; is the lift coefficient; and Cp 
is the drag coefficient. Note that, for the aeroassisted flight experiment, 
the angle of attack is kept constant; the aerodynamic coefficients are assumed 
to be independent of the Mach number and the Reynolds number; and the spacecraft 
is controlled via the angle of bank. 
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9. Transformation Relations 

In this section, we supply some transformation relations which allow 
one to pass from (i) quantities computed in an Earth-fixed system to 
(ii) quantities computed in an inertial system, and viceversa. 

9.1. Spacecraft Position . Let r, 9, ({> denote the spherical coordinates of 
the spacecraft P in the Earth-fixed system Ox e y g z e . Let r^ , 9.., <J>.. 
denote the spherical coordinates of the same spacecraft in the inertial 
system Ox.y.z^. Assume that the axes of the Earth-fixed system coincide 
with the axes of the inertial system at time instant t = 0. Then, the 
following transformation relations hold: 


r i = r ’ 

(66a) 

0,. = 0 + U)t, 

(66b) 

$.j = <f>- 

(66c) 

Equations (66) imply the following inverse relations: 
r = r., 

(67a) 

9=9^ - wt. 

(67b) 

<j) = (f)^ . 

(67c) 

9.2. Spacecraft Velocity. Let V, y, x denote the 

velocity modulus, the 


path inclination, and the heading angle in the Earth-fixed system 0x e y e z g . 
Let V., y . , x-j denote the velocity modulus, the path inclination, and the 
heading angle in the inertial system Ox^y^z-. Let V denote the velocity 
vector in the Earth-fixed system; and let V. denote the velocity vector 
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in the inertial system. 

We employ the theorem of composition of velocities, which states 

-y 

that the inertial velocity is the sum of the relative velocity V 
and the transport velocity to x OP, 

-y -y *+* 

V i = V + to x OP. (68) 


The vectors appearing in Eq. (68) can be written in terms of their 
components on the local horizon system p *h y h z h as follows [see Eqs. (48), 
(50), (51)]: 

->■ -*■ -*• 

V i = V^osy. cosx i i h + V i cosy i sinx i J h - V^iny. k h> (69a) 

-y -*■ ■+ 

V = Vcosy cosy i^ + Vcosy sinx - Vsiny k^, (69b) 

coxOP = corcos<}> i h> (69c) 


with the implication that 

V^cosy^ cosx^ = Vcosy cosy + tor cos<j>, (70a) 

V^osy^ sinx^ = Vcosy siny, (70b) 

V.siny. = Vsiny. (70c) 

Laborious manipulations, omitted for the sake of brevity, lead to the 
following transformation relations: 

A 2 

V.j = /[v + ZtorVcosy cosy cos<j> + (cor cos<}>) ], 


(71a) 
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2 2 

tany^ = Vsiny/yUVcosy) + 2wrVcosy cosx cos<f> + (wr cos4>) ], (71b) 

tanx.j = V cosy sinx/(Vcosy cosy + wrcos<j>). (71c) 

Equations (71) imply the following inverse relations: 

2 2 

V = /[V.. - 2wr.jV..cosY.j cosx i cost^^. + (wr .coscj)^ ) ], (72a) 

2 2 

tany = V 1 .siny i //[(V 1 .cosy^ ) - 2wr i -V^cosy i cosy.,- cos^^ + (w^cost^. ) ], (72b) 

tanx = V^cosy^ sinx i /(V i cosy i cosx^ - wr..cos4>.j ) . (72c) 

9.3. Cartesian Coordinates . After the spacecraft position is known in spherical 
coordinates, the corresponding Cartesian coordinates can be computed. The 
following transformation relations hold: 

x g = rcose cos<J> = r^cos(9 i - wt)cos<j>.., (73a) 

y g = rsine cos<j> = r^sin(0^ - wt)cos<t>.j, (73b) 

z g = rsi n4> = r^sin^^,. (73c) 

Equations (73) imply the following inverse relations: 

x. = r^cose^ cos^ = rcos(6 + wt)cos<(), (74a) 

y. . = r^sin©^ cosc^. = rsin(0 + wt)cos<f>, (74b) 


= rsint}). 


(74c) 
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10. Conclusions 

This report is the first of a series dealing with the determination 
of optimal trajectories for the aeroassisted flight experiment (AFE). The 
AFE refers to the study of the free flight of an autonomous spacecraft, 
shuttle-launched and shuttle-recovered. Its purpose is to gather atmospheric 
entry environmental data for use in designing aeroassisted orbital transfer 
vehicles (AOTV). 

It is assumed that: the spacecraft is a particle of constant mass; 
the Earth is rotating with constant angular velocity; the Earth is an oblate 
planet, and the gravitational potential depends on both the radial distance 
and the latitude; however, harmonics of order higher than four are ignored; 
the atmosphere is at rest with respect to the Earth. 

Under the above assumptions, the equations of motion for hypervelocity 
atmospheric flight (which can be used not only for AFE problems, but also 
for AOT problems and space shuttle problems) are derived in an Earth-fixed 
system. Transformation relations are supplied which allow one to pass from 
quantities computed in an Earth-fixed system to quantities computed in an 
inertial system, and viceversa. 
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